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a: 

q I Abstract. Motivated by the general problem of studying sample-to-sample 

'— '■ fluctuations in disorder-generated multifractal patterns we attempt to investigate 

I . analytically as well as numerically the statistics of high values of the simplest model 

^ 1 - the ideal periodic 1/f Gaussian noise. Our main object of interest is the number 

of points Mm(x) above a level §U m , with V m = 2 InM standing for the leading-order 
typical value of the absolute maximum for the sample of M points. By employing the 
■ thermodynamic formalism we predict the characteristic scale and the precise scaling 

. form of the distribution of Nm{x) for < x < 2. We demonstrate that the powerlaw 

" forward tail of the probability density, with exponent controlled by the level x, results 

in an important difference between the mean and the typical values of Mm{x). This 
' . I , can be further used to determine the typical threshold x rn of extreme values in the 

pattern which turns out to be given by im P ' = 2 — clnln M/bxM with c = |. 
. Such observation provides a rather compelling explanation of the mechanism behind 

$H I universality of c. Revealed mechanisms are conjectured to retain their qualitative 

validity for a broad class of disorder-generated multifractal fields. In particular, 
we predict that the typical value of the maximum p ma x of intensity is to be given 
by —\np max = a_ In A/ + 2f(a J mm -^ + 0(1), where f(a) is the corresponding 
singularity spectrum vanishing at a = a- > 0. For the 1/f noise case we further study 
asymptotic values of the prefactors in scaling laws for the moments of the counting 
function. Our numerics shows however that one needs prohibitively large sample sizes 
to reach such asymptotics even with a moderate precision. This motivates us to derive 
exact as well as well-controlled approximate formulas for the mean and the variance 
of the counting function without recourse to the thermodynamic formalism. 
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1. Introduction 

Investigations of multifractal structures of diverse origin is for several decades a very 
active field of research in various branches of applied mathematical sciences like chaos 
theory, geophysics and oceanology [H [2] as well as climate studies [3], mathematical 
finance [HE], and in such areas of physics as turbulence [E1I7], growth processes [8], and 
theory of quantum disordered systems [9]. The main characteristics of a multifractal 
pattern of data is to possess high variability over a wide range of space or time scales, 
associated with huge fluctuations in intensity which can be visually detected. Another 
common feature is presence of certain long-ranged powerlaw-type correlations in data 
values [TO]. 

To set the notations, consider a certain (e.g. hypercubic) lattice of linear extent 
L and lattice spacing a in d— dimensional space, with M ~ (L/a) d ^> 1 standing for 
the total number of sites in the lattice. The multifractal patterns are then usually 
associated with a set of non- negative "heights" hi > attributed to every lattice site 
i — 1, 2, . . . , M such that the heights scale in the limit M — > oo differently at different 
sites: hi ~ M x% [§[ with exponents Xi forming a dense set. To characterize such a 
pattern of heights quantitatively it is natural to count the sites with the same scaling 
behaviour. Then a multifractal measure is characterized by a (usually, concave) single- 
maximum singularity spectrum function f(x). Denoting the position of its maximum 
as x = Xq } such function describes the (large-deviation) scaling of the number of points 
in the pattern whose local exponents x% belong to some interval around xq- More 
precisely, defining the density of exponents by Pm{x) = J2iLi $ (nTM — x ) a nontrivial 
multifractality implies that such density should behave in the large-M limit as [11] 



with a prefactor cm(%) of the order of unity which may still depend on x. We will refer 
below to the above form as the multifractal ansatz. The major effort in the last decades 
was directed towards determining the shape and properties of f(x). In contrast, our 
main object of interest will be the behaviour of the prefactor cm{x) which is much 
less studied, to the best of our knowledge. In particular, if the multifractal pattern is 
randomly generated like e.g. those considered in [9] the ansatz ([[]) is expected to be 
valid in every realization of the disorder. One may then be interested in understanding 
the sample-to-sample fluctuations of the prefactor cm(x). 

§ Usually one defines exponents on via the relation hi ~ L ai i.e. by the reference to linear scale 
L instead of the total number of sites M ~ (L/a) d , and similarly for the density of exponents 
p(a) ~ L^ a \ We however find it more convenient to use instead the exponents Xi = ai/d and 
the singularity spectrum f(x) — \f{a)- 
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Further introducing the counting functions 

POO 

N > (x)= p M (y)dy, N < (x)= I PM (y)dy (2) 



for the total number N > (x) of sites of the lattice where heights satisfy hi > M x 
(respectively, hi < M x ), substituting the multifractal form of the density into (j2J) and 
performing at InM ^> 1 the resulting integral for x > Xq by the Laplace method we 
find N > (x) « cm (x)M^ x ^ / \ f (x) | \/ln M and a similar expression for N<(x) for x < xq, 
relating the singularity spectrum f(x) to the counting functions. As both N > (x) and 
N<(x) can not be smaller than unity we necessarily have f(x) > for all x, and the 
condition f(x) = defines generically the maximal x + and the minimal x_ threshold 
values of the exponents which can be observed in a given height pattern. 

The singularity spectrum f(x) is not a quantity which is easily calculated 
analytically or even numerically for a given multifractal pattern of heightspj]. An 
alternative procedure of analysing the multifractality is frequently referred to in the 
literature as the thermodynamic formalism [H12]- In that approach one characterizes the 
multifractal pattern by the set of exponents ( q describing the large-M scaling behaviour 
of the so-called partition functions Z q as 

M 

Z q = J2 h !~ M<q > lnM>l (3) 

i=i 

To relate ( q to the singularity spectrum f(x) discussed above one rewrites ([3]) in terms 
of the density as Z q = J_ M qx Pm{x) dx, and again employs the multifractal ansatz 
(JIJ) for pm(x). Evaluating the integral in the InM > 1 limit by the steepest descent 
(Laplace) method gives 

( 2tt \ 1/2 

z q ~ c M (x*) \ TjF(Z )| J MCq wliere f'( x *) = -q and Cq = f{%*) + qx*, (4) 

where we have assumed that x_ < < x + for simplicity. This shows that the relation 
between ( q and f(x) is given essentially by the Legendre transform. We thus see that 
formally the original definition of multifractality based on the density (or, equivalently, 
the counting functions N >j< (x)) and the thermodynamic formalism approach ()3])- (J4j) 
should have exactly the same content for InM — > oo, provided the singularity spectrum 
is concaveJJJ]. Note also the normalization identity Z = Pm(u) dy = M implying 
Co = 1- It also shows that at the point of maximum x = Xq we must necessarily have 
f(xo) = 1 and that cm(xq) is indeed of the order of unity. 

In practice, to extract singularity spectra from a given multifractal pattern 
obtained in real or computer experiments, one frequently employs the so-called box 

| Examples of non-concave multifractality spectrum and the associated thermodynamic formalism 
are discussed in [12] 
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counting procedure which can be briefly described as follows. Subdivide the sample 
into Mi = (L/l) d non-overlapping hypercubic boxes of linear dimension I. Associate 
with each box the mean height Hk(l) = {l/a)~ d J2ien k hi, and define the scale-dependent 
partition functions 

^,i) = f E^f ( 5 ) 

1 k=l 

Note that for / = a obviously M a = M = (L/a) d and Z q (a,L) coincides with the 
partition function Z q featuring in the thermodynamic formalism. One may however 
observe that in the range a <^ I <^ L the scale-dependent partition functions are 
sensitive to the spatial correlations in the heights at different lattice points. In 
particular, a simple consideration shows that when the heights are powerlaw-correlated 
in space as is actually the case for the most of systems of interest [10] (see also below) 
the scaling behaviour of Z q (l,L) depends non-trivially on both I /a and L/a. At the 
same time the behaviour of the combination I g (l, L) = Z q {l, V)j [Zi(l, L)] q turns out to 
be a function only on the scaling ratio L/l and is given by 

IS, I) = ~ (f ) ' ' where r * = d ^ - Cr) ( 6 ) 

which allows to get reliable numerical values of the scaling exponents r q by varying 
the ratio L/l over a big range. Further noticing that for q = 1 the /—dependence of 
the partition function disappears due to linearity: Zi(l,L) = (1/ a)~ d (L/l)~ d ~Y/^ hi ~ 
(L/a)^~ d we also can reliably extract Ci from the same data, hence relate the set of 
exponents r q to ( q for q ^ 1. 

The quantities I q = I q (a, L) have interpretation of the inverse participation ratio's 
(IPR's) and are very popular in the theory of the Anderson localization [9] and related 
studies. Passing from the partition functions Z q of the thermodynamic formalism to the 
IPR's is equivalent to focusing on the properties of the normalized probability measure 
< pi = hi/Zi < 1, Pi = 1 rather than on the original height pattern hi itself. In fact 
in such a setting it is more natural to introduce the scaling of those weights in the form 
Pi ~ M~ a \ on > and consider the corresponding singularity spectrum f(a) related 
directly to the Legendre transform of the exponents r q . Working with the exponents r q 
has some advantages, as one can show they must be monotonically increasing convex 
function of q: ^ > 0, ^ < 0. In many situations, as e.g. the diffusion-limited 
aggregation [5] or indeed the Anderson localization the multifractal probability measures 
arise very naturally. In other contexts, e.g. in turbulence or in financial data analysis, 
the normalization condition seems superfluous. In the main part of the present paper 
we are mainly interested in the pattern of heights and are therefore concentrating on 
partition functions. We will discuss the normalized multifractal probability measures 
in the end. 
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The above description is generally valid for a given single multifractal pattern of 
any nature, not necessarily random. In recent years considerable efforts were directed 
towards understanding disorder-generated multifractality, see e.g. [9j [131 ttH 02] for 
a comprehensive discussion in the context of Anderson localisation transitions and 
various associated random matrix models, [16], [17] in the context of harmonic measure 
generated by conformally invariant two-dimensional random curves and [HJ [191 EQ] 
for examples related to Statistical Mechanics in disordered media. We just briefly 
mention here that one of the specific features of multifractality in the presence of 
disorder is a possibility of existence of two sets of exponents, r q versus f g , governing 
the scaling behaviour of the typical IPR denoted Iq ~ M~ Tq versus disorder averaged 
("annealed") IPR, I q ~ M~ Tq . Here and henceforth the overline stands for the averaging 
over different realisations of the disorder. Namely, it was found that for large enough 
q > q c the two exponents will have in general different values: r q ^ f q . The possibility of 
"annealed" average to produce results different from typical is related to a possibility of 
disorder-averaged moments to be dominated by exponentially rare configurations. As a 
result, the part of the "annealed" multifractality spectrum recovered via the Legendre 
transform from f q for q > q c will be negative [21], [9]: f(x) < for x < x_, and similarly 
for x > x + . Further detail can be found in the cited papers and in the lectures [22J. 

Another important aspect of random multifractals revealed originally by Mirlin 
and Evers [14] in the context of the Anderson localization transition is the fact that 
IPR's I q for disorder- induced multifractal probability measures are generically power- 
law distributed: V(I q /I { q t] ) ~ {I q / if 1 )' 1 -"* [H ^. Such behaviour suggests that 
the actual values of the counting functions iV> (x) , N< (x) should also show substantial 
sample-to-sample fluctuations, even in the range x_ < x < x + where the singularity 
spectrum f(x) is self- averaging and the same multifractal scalings N > (x) ~ M^ x ' is to 
be observed in every realization of the pattern. Though the presence of such fluctuations 
was already mentioned in [18], a detailed quantitative analysis seems to be not available 
yet. The main goal of our paper is to achieve a better understanding of statistics of 
the counting functions N > (x), N < (x) by performing a detailed analytical as well as 
numerical study of arguably the simplest, yet important class of multifractal disordered 
patterns - those generated by one-dimensional Gaussian processes with logarithmic 
correlations, the so-called 1/f noises. 

The structure of the paper is as follows. In the next section we will introduce 
the 1/f noise signals and discuss their properties already known from the previous 
works. Then we will use that knowledge to show that the probability density of the 
counting function N > (x) for such a model is characterized by a limiting scaling law with 
a powerlaw forwards tail, with the power governing the decay changing with the level 
x. We will then demonstrate that such powerlaw decay has nontrivial implications for 
the position of the maxima (or, with due modifications, minima) of such processes, and 
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derive the expression for the threshold of extreme values. Finally, using 1/ / noises as 
guiding example we will attempt to reinterpret the results of the theory developed in 
[14J to get a rather general prediction for the position of extreme value threshold for a 
broad class of disorder-generated multifractal patterns whose intensity is characterized 
by power-law correlations. We conclude with briefly discussing a few open questions. 

2. 1/f noise: mathematical model and previous results 

An ideal 1/f (or "pink") noise is a random signal such that spectral power (defined 
via the Fourier transform of the autocorrelation function of the signal) associated 
with a given Fourier harmonic is inversely proportional to the frequency u = 2-hf. 
Signals of similar sort are known for about eighty years and believed to be ubiquitous 
in Nature, see [23] for a discussion and further references. Rather accurate 1/f 
dependences may extend for several decades in frequency in some instances, es. e.g. 
in voltage fluctuations in thin-film resistors [23] or resistance fluctuations in single- 
layer graphene films [25], in non-equilibrium phase transitions [26] . and in spontaneous 
brain activity |27j. Still, the physical mechanisms behind such a behaviour are not 
yet fully known, and are a matter of active research and debate. It was noticed quite 
long ago (see e.g. [28]) that a generic feature of all such signals is that the two- 
point correlations (covariances) depend logarithmically on the time separation. During 
the last decade it became clear that random functions of such type appear in many 
interesting problems of quite different nature, featuring in physics of disordered systems 
[18j,[29j,[30],[3T],[32], quantum chaos [33], mathematical finance [331 EH], turbulence 
[361 EH EE] and related models [391 SOI 131] , as well as in mathematical studies of random 
conformal curves [32], Gaussian Free Field [331 [33] and related models inspired by 
applications is statistical mechanics [32] and quantum gravity [46] . and the most recently 
in the value distribution of the characteristic polynomials of random matrices and the 
Riemann zeta- function along the critical axis [47] . Let us note that a simple argument 
outlined in [22J and repeated in section 5 of the present paper shows that by taking the 
logarithm of any spatially homogeneous powerlaw-correlated multifractal random field 
we necessarily obtain a field logarithmically correlated in space. A somewhat similar 
in spirit suggestion to refocus the attention from the multifractal ("intermittent") 
signals and fields to their logarithms was also put forward in [3?]. All this makes 
logarithmically-correlated Gaussian processes an ideal laboratory for studying disorder- 
induced multifractality, though investigating the effects of non-Gaussianity remains a 
challenging outstanding issue. 

Despite the fact of being of intrinsic interest and fundamental importance, a 
coherent and comprehensive description of statistical characteristics of ideal 1/f noises 
seems not to be yet available, and relatively few properties are firmly established even for 
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the simplest case of a Gaussian 1/f noise. Among the works which deserve mentioning 
in such a context is the paper [48] which provided an explicit distribution of the "width" 
(or "roughness") for such a signal, as well as the work [49] describing a curious property 
of spectral invariance with respect to amplitude truncation. In recent papers [30].[3T] 
and [45] the statistics of the extreme (minimal or maximal) values of various versions 
of the ideal 1/f signals was thoroughly addressed. From that angle the subject of the 
present paper is to provide a fairly detailed picture of statistics of the number of points 
in such signals which lie above a given threshold set at some high value. The latter can 
be rather naturally defined as being at finite ratio to the typical value of the absolute 
maximum. 

In this paper we are going to consider only Gaussian ideal 1/f noises, 27r-periodic 
version of which is naturally defined via a random Fourier series of the form 



where v n is a set of i.i.d. complex Gaussian variables with zero mean and the variance 
\v n \ 2 = 1, with the asterisk standing for the complex conjugation, and the bar for the 
statistical averaging. It implies the following covariance structure of the random signal 



Mathematically such series represents the periodic version of the fractional Brownian 
motion with the Hurst index H = 0. The corresponding definition is formal, as the series 
in (J7|) does not converge pointwise, the fact reflected, in particular, in the logarithmic 
divergence of the covariance in (JEW Although it is possible to provide several bona 
fide mathematically correct definitions of the ideal 1/f noise as a random generalized 
function (based, for example, on sampling 2c? Gaussian free field along the specified 
curves , e.g. the unit circle for the periodic noise, see [42], or the constructions proposed 
in [39] or [45]), for all practical purposes the 1/f noises should be understood after a 
proper regularization. In what follows we will use explicitly the regularization proposed 
by Fyodorov and Bouchaud [30] , though we expect the main results must hold, mutatis 
mutandis, for any other regularization. 

In the model proposed in [30] one subdivides the interval t G (0, 2tt] by finite 
number M of observation points = jj-k where fc = l,...,M<oo, and replaces the 
function V(t), t G [0, 2ir) with a sequence of M random mean-zero Gaussian variables 

% One can also define other, in general non-periodic versions of similar log-correlated random processes 
on finite intervals using different basis of orthogonal functions, or even exploit the appropriate random 
Fourier integral to define the process on the half-line < t < oo. The corresponding models arise very 
naturally in the context of Random Matrix Theory and will be discussed in a separate publication 50 




n=l 



(7) 
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Vk correlated according to the M x M covariance matrix Ck m — VkV m such that the 
off-diagonal entries are given by 

C k ^ m = -21n\2sm^{k-m)\. (9) 

To have a well-defined set of the Gaussian-distributed random variables one has to 
ensure the positive definiteness of the covariance matrix by choosing the appropriate 
diagonal entries Ckk- A simple calculation (30] shows that as long as we choose 

C kk = V% >21nM, Vfc = l,...,M. (10) 

the model is well defined, and we will actually take the minimal possible value: 
Ckk = 21nM + e, Wk with a small positive e < 1. We expect that the statistical 
properties of the sequence Vj. generated in this way reflect correctly the universal 
features of the 1/f noise. An example of the signal generated for M = 4096 according 
to the prescription above via the Fast Fourier Transform (FFT) method as explained 
in detail in [31] is given in the figure. 




Figure 1. A sample of the regularized periodic 1/f noise for M = 4096 observation 
points. The upper broken line marks the typical value of the maximum V m = 
2 In M — | In In M and black dots mark a few exceedances of that level. The lower 
line is at the level -^Vm and the blue dots at the bottom mark points i supporting 

Vi > -7=V m . The set looks like a fractal. 

v 2 

Using such a model, the authors of [30] defined the associated random energy model 
via the partition function Z{f3) = e_/?V N with the temperature T = > 

and succeeded in determining the distribution of Z{j5) in the range < 1. To 
reinterpret those findings in the context of multifractality we introduce the height 
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variables hi = e Vi > and rename (3 — > —q , converting Z[f3) of the random energy 
model to Z q of the "thermodynamic formalism", eq.(j3]). Note that due to the statistical 
equivalence of Vi and — Vi in the model all results may depend only on \q\. Then the 
findings of [30J can be summarized as follows. The probability density of the random 
variable Z\ q \<i consists of two pieces, the body and the far tail. The body of the 
distribution has a pronounced maximum at Z ~ Z e (q) = M 1+q /Y(l - q 2 ) < M 2 , 
and a powerlaw decay when Z e <Z< M 2 . Introducing z = Z q /Z e (q) the probability 
density of such variable is given explicitly by 

i+ 1 i 

Vq{z) = Y 2 iS) ^ e_( * P ' ^« Ml "" 2 ' M <l ( u ) 

For z 3> M l ~ q2 the above expression is replaced by a lognormal tail|30j. Note that the 
probability density ( TTTj) is characterized by the moments 

?m»i ~ T(l - sq 2 ), Re s < q- 2 (12) 
where T(z) is the Euler Gamma-function. It is worth noting that although the particular 
form of the density ( flTI) is specific for the chosen model of 1/f noise, the power-law 
forward tail V q {Z) oc Z c ' 2 is expected to be universal (31] and so the divergence 
of moments of the partition function for Res > q~ 2 . A closely related fact which can 
be also traced back to the existence of the universal forward tail is that the typical 
scale Z e (q) behaves for q — > 1 as Z e (q) /M 1+q2 ~ (1 — q) — > 0. This property will have 
important consequences at the level of the counting function. 

To get some understanding of how the above asymptotic results agree with the 
direct numerical simulations of the model for large, but finite M it is useful to provide 
the exact finite-M expression for the second moment of the partition function [30] 

Z 2 = M 1+4q2 + M 2 ^ Si M \q% SfV) = ^ E H* 2 (£0r (13) 

i=i 

Second term here is dominant in the large-M limit and gives precisely the result ( TT2l) 
as asymptotically we can replace the sum by the integral convergent for q 2 < 1/2 and 
get 

lim S[ M \q 2 ) = 2 - f ^ [4sin^]- 2 d9 = (14) 

The main correction to this asymptotic result is given by the first term in (Tl3|) . whose 
relative contribution is small as M l ~ 2q2 for q 2 < 1/2. In figure |2] we show the exact 
evaluation for the relative variance 5 z (q, M) = -0^ — 1 with z = Z q /Z e (q) given explicitly 
by 

5 z (q,M) = S { 2 M \q 2 ) + M 2q2 - 1 -l (15) 

We see that 5 z (q, M) clearly approaches the asymptotic value 5 z (q, oo) = [w^zpsA — 1 
for M ~ 2 11 . 
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Figure 2. Convergence of the relative variance to the asymptotic value. Red dots 
correspond to q = 1/3 and blue dots correspond to q = 1/2. The asymptotic behavior 
is reached for M ~ 2 11 . When q increases the convergence slows down. 



3. Thermodynamic formalism for the counting function of the 1/f noise 
sequence and the threshold of extreme values 

The statistics of Z g for the model under consideration suggest that it reflects the 
corresponding strong sample-to-sample fluctuations in the counting function of pattern 
of heights. Our goal is to quantify statistics of those fluctuations by considering the 
total number i\T>(a;), which in the present context will be denoted as Nm{x), of the 
x—high points in the (regularized) 1/f sequence Vi, . . . , Vm such that Vi > x In M which 
is equivalent to hi = e Vl > M x . 

We then relate the number Nm{%) to the partition function Z q by the 
thermodynamic formalism: 

poo /»oo 

Mm{x) = In M / p M (y) dy, Z q = \nM M™p M (y) dy (16) 

J X J — oo 

where now the density Pm{v) = = S/tli ^(^fc ~ U^ n M) is anticipated to be given 
in the large-M limit by the multifractal ansatz of the particular "improved" form 

PM(y) * c M (y)-j=-, c ^) = 20Fr( l-,V4) ; M<2 (17) 
with a random coefficient nuiv) of order of unity which strongly fluctuates from 
one realization of the sequence Vi to the other in such a way that its probability 
density is given by (ITT]) with q = y/2. Indeed, substituting such form of the density 
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to Z q in ( 1161) and performing the integrals in the limit InM ^> 1 by the Laplace 
method we arrive at the asymptotic behaviour Z q « n M {2q) Z e (q), q < 1, where 

Z e {(l) = r(i-q 2 ) P rec i se ly as we have found from the exact solution ffTTj) . On the other 
hand, substituting the same ansatz to the counting function in ( ll~6l) yields by the same 
method A/m>i(x) ~ n>M( x ) Nt(x) where the typical value Nt(x) is given by 

M l-* a /4 J 



X ) 



< x < 2. 



xVvrlnMr(l -x 2 /4) : 
Thus, our main conclusion is that the random variables n = Mm{x) / Mt{x) and 
z = Z q /2/Z e (q/2) must be distributed in the large— M limit according to the same 
probability law, which after invoking (TTTj) yields the asymptotic probability density for 
the scaled counting function in the form 

V x (n) = -^-^-r f 0<x<2. (19) 



r 2 i+ 



4 

7? 



The shape of the distribution for a few values of x is presented in Fig. [31 



c 

x 
Q_ 




Figure 3. From top to bottom the probability density Eq.(|19p for x = 1, x = 1.4, 
a; = 1.7. Arrows indicates the position of the mean, the typical value corresponds 
always to 1 in the chosen scaling. 



The following qualification is needed here. For large but finite M such form of the 
density stops to hold true for extremely large n — > oo as in any realization obviously 
Nm(x) < M at the very least. Therefore there must exist an upper cut-off value Af c (x) 
such that for n > n c = Af c (x) / J\ft(x) the scaling form of the probability density f lT9|) 
loses its validity. Similarly, another restriction on the validity of ( !T9l) should exist in the 
region of extremely small n — > due to implicit condition Nm{%) ^ 1- Precise value 
of the cutoffs can not be extracted in the framework of the thermodynamic formalism, 
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but in the end of the paper we will get an estimate for n c in the range a/2 < x < 2 
from independent calculations, see f )32|) . 

The important scale J\f t (x) defined explicitly in ( fl8l) describes typical values of the 
counting function J\f M (x) for a given observation level x. In particular, it can be used 
to define one of the objects of central interest in the present paper, the threshold of 
extreme values. The latter stands for such a level above which typically we can find for 
In M ^> 1 only a few, i.e. of the order of one points of our random sequence. The scaling 
behaviour Nt(x) ~ M^ x \ f(x) = 1 — x 2 /4 is the hallmark of the multifractality. A 
very similar parabolic singularity spectrum characterizes the high value pattern of the 
two-dimensional Gaussian free field as revealed in [18] and proved in a mathematically 
rigorous way in [H]. The above result for f(x) is the simple one-dimensional analogue 
of that fact, see e.g. [15] where it is rigorously shown that Mitim^^ ln ^^ = 1 — x 2 /A 
in our notations. Note that as the singularity spectrum f(x) vanishes at x = 2 the 
typical position of the absolute maximum of the random sequence of ViS is given by 
V max = 2 In M at the leading order. The corresponding subleading term was conjectured 
in the work by Carpentier and Le Doussal [29] to be V m = 21nM — clnlnM + 0(1), 
with c = 3/2. That conclusion was based upon an analysis of the travelling wave-type 
equation [51] appearing in the course of one-loop renormalization group calculation, 
and the value c = 3/2 was conjectured to be universally shared by all systems with 
logarithmic correlation. Such a result is markedly different from c = 1/2 typical for 
short-ranged correlated random signals, so the value of c may be used as a sensitive 
indicator of the universality class. Indeed, in a recent numerical studies of the behaviour 
of the logarithm of the modulus of the Riemann zeta-function along the critical line 
[17] the value 3/2 was used to confirm the consistency of describing that function as 
a representative of logarithmically correlated processes. Despite its importance, no 
transparent qualitative argument explaining c = 3/2 vs. c — 1/2 values was ever 
provided, to the best of our knowledge, though for the case of the 2D Gaussian free 
field the value 3/2 was very recently rigorously proved by Bramson and Zeitouni [52J 
by exploiting elaborate probabilistic arguments. Below we suggest a very general and 
transparent argument showing that the change from c = 1/2 to c = 3/2 is a direct 
consequence of the strong fluctuations in the counting function reflected in the power- 
law decay of the probability density ( 1191) . That observation not only allows one to 
explain the origin of c = 3/2 for the Gaussian case, but has also a predictive power in 
a more general situation as will be demonstrated in the section 5. 

To begin with presenting the essence of our argument we first observe that ffT9]) 
implies that n = T ^1 — ^\ , < x < 2 so that the mean value of the counting function 
is given asymptotically by 
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We shall see in the next section that the above expression is asymptotically exact for any 
real x > without restriction to < x < 2. Notice however that the mean value ( 120]) 
and the characteristic scale N t (x) differ from each other by the factor T (i-x 2 /4) tending 
to zero as x — > 2. Such a difference, origin of which can be again traced to the specific 
power-law tail of the probability density (Tl9|) . see Fig. [3] and is one of the hallmarks of 
the random signals with logarithmic correlations. Indeed, consider for comparison the 
case of uncorrelated i.i.d. Gaussian sequence sharing the same variance Vf = 21nM 
with the logarithmically correlated noise (by historical reasons it is natural to refer to 
such model as the Random Energy Model, or REM [53]). A straightforward calculation 
shows that we still would have precisely the same mean value of the counting function 
as in the logarithmic model, but unlike the latter it will be simultaneously the typical 
value of that random variable as no powerlaw tail is present in that case (see Fig. H] 
and discussion in the next section). 

Such a difference between the two cases has important implications for the location 
of the threshold x = x m which corresponds to the region of extreme value statistics of 
multifractal heights. Indeed, by approximating the singularity spectrum close to its 
right zero x + = 2 as f(x) ~ (2 — x), and similarly writing T ^1 — oc (2 — x)^ 1 we 
observe that the condition J\f t (x m ) ~ 1 is equivalent to the equation: 

(2 -x m )\nM- -lnlnM + ln(2 - x m ) = (21) 

solving which for In M ^> 1 to the first non-trivial order gives precisely x m = 2 — c lr \^^f 
with c = 3/2. Had we replaced in the above condition the typical value Af t (x m ) with 
the mean A/Af(x m ) we would arrive to the same expression for x m but with the value 
c = 1/2 replacing c = 3/2. This perfectly agrees with such x m being the extreme value 
threshold for short-ranged correlated random sequences. The suggested explanation of 
the transmutation of the coefficient c based on the "typical versus mean" argument 
seems to us very transparent and supports the conjectured universality of the result. 
Indeed, the thermodynamic formalism combined with the results of [3T] suggests that 

4 

the power- law forward tail V x (Af) oc i (77-)^ should be universal for one-dimensional 
Gaussian processes with logarithmic correlations. It also should show up, mutatis 
mutandis, in higher- dimensional versions of the model like e.g. the Gaussian Free 
Field on the lattice. Such a tail will ensure the difference between the typical and the 
mean of the counting function by a factor which vanishes linearly on approaching the 
extreme value threshold value x = x m . This factor then will lead to the value c = 3/2 by 
the mechanism illustrated above for our particular explicit example. The case of non- 
Gaussian signals with logarithmic correlations is relevant for general disorder-generated 
multifractals and is discussed in Section 5 of the paper. 

Finally it is worth mentioning that the fact that for x > 2 the mean value of the 
counting function (|20|) is exponentially small reflects the need to generate exponentially 
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large number of samples to have at least a single event with Vi > 2x In M for In M ^> 1 
as such values of Vi will not show up in a typical realization (cf. earlier discussion about 
"annealed" vs. "quenched" singularity spectra). 

4. Exact results versus asymptotics 

The above results for the counting function obtained in the framework of the 
thermodynamic formalism are expected to be valid as long as InM ^> 1. To get a 
feeling of how big in practice should be M to ensure the validity of our asymptotic 
formulae, it is natural again to try to perform the direct numerical simulations of 
the regularized version of the ideal 1/f noise. We start with checking directly the 
distribution of the scaled counting function, Eq.( Tl9]) . for a particular values x — 1. 
The results are presented in FigJH They show that although the main qualitative 
features of the distribution (in particular, the well-developed powerlaw tail) are clearly 
in agreement with the theoretical predictions, the curve is still rather far from its 
predicted asymptotic shape for M = 2 18 , and the convergence is too slow to claim a 
quantitative agreement. 




0.5 1 1.5 2 2.5 3 3.5 4 



n 

Figure 4. Probability density for the scaled counting function, V x (n) with x = 1. 
Black solid line corresponds to the analytical prediction EqOlJ] Red lines correspond 
to numerical simulation of the regularized 1/f Gaussian signal generated via the 
Fast Fourier Transform (FFT) method (lower red curve M — 2 14 , higher M — 2 18 , 
data collected from 10 5 samples.). The blue line corresponds to the REM model: 
M = 2 18 i.i.d. zero mean Gaussian variables with variance 21nM; the density for 
n = Nm{x)/Nm(x) is clearly converging to the delta peak and does not show any 
power-law tails, see discussion in the text. 



To get a better understanding of the mechanism of such disagreement at a 
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quantitative level, and to check the results obtained in the framework of the 
thermodynamic formalism we choose to consider in much greater detail the first two 
moments of the counting function. The asymptotic formula ( 1T91) yields for the mean of 
the counting function the expression (1201) and for the variance 

-1 2 

r(i-x 2 /2) 



Wm{x)\ 



.-, 2 



r 2 (i-x 2 /4) 



- 1, < x < V2 



(22) 



At the same time it is possible to derive a closed form exact finite M expression for 
the first two moments of the counting function Nm{x) = Y^tLi Q(Vi ~ xlnM) without 
any recourse to the thermodynamic formalism. Here we have used the Heaviside step 
function 6{u) = 1 for u > and 9{u) = otherwise. The mean value can be immediately 
computed as 



M 



>x\nM 



cxp 



V" 



41nM 



alv 



M 



Erfc (-VlnM 



(23) 



and is independent of the correlations. The problem of deriving a closed-form expression 
for the variance which is amenable to accurate numerical evaluation for very big 
InM ^> 1 is less trivial and may have an independent interest. Before presenting 
the results we find it most convenient to define the following object 



A M (x) = N M (x) (N M (x) - 1) 



M- 1 

M 



n 2 



Nm(x) 



in terms of which the relative variance is expressed as 

-1 2 

1 1 

_ M 4 



[Hm{x)Y- M m {x) 



N m (x) 



n 2 



M M {.x) 



+ 5 n (x;M), 5 n {x;M) 



^m{x) 



N M {.x) 



(24) 



(25) 



Am{x) is a convenient measure of correlation- induced fluctuations. Using the definition 
of the counting function we explicitly get: 

A M (x) = [W - x ln MfflYj ~ x In M) - 6(Vi - x In M) ■ 9(Vj — x In M) (26) 

from which is evident that Am{x) vanishes for any i.i.d. sequence. In the latter case the 
relative variance tends to zero in the large-M limit assuming Mm{x) — > oo for M — > oo. 
This simply means that in the i.i..d. case the variable Mm{x) / 'Nm{x) is self-averaging, 
i.e its limiting density approaches the Dirac delta-function, see Fig. |H On the other 
hand, Am{x) is formally different from zero for correlated variables. Nevertheless, 
using the general formalism exposed in the Appendix A one can satisfy oneself that 
for stationary Gaussian sequences with correlations decaying fast enough to zero at big 
separations (e.g. as a power of the distance) the quantity 5 n (x; M) still tends to zero 
as M — > oo. In contrast, we will see below that in the logarithmically correlated case 
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5 n (x; M) tends to a finite positive number for M — > oo, and thus coincides with the 
leading behaviour of the variance of the counting function. 

In the Appendix A we have derived the exact expression for Am(x) for a general 
correlated Gaussian sequence. For the periodic 1/ / noise sequence using (JSJ the result 
reads 



In the figure we test its validity by comparing the results obtained for moderate 
value of M by direct numerical simulations of the log-correlated sequences with the 
predictions of (l2Tj) . 

The expressions above are suitable for developing a well-controlled approximation 
to the exact expression ( 12 7|) in the large-M limit assuming lnM> 1. First of all, it is 
clear that in the large-M limit we may replace the discrete sum in fl27|) by the integral 
treating as a continuous variable 9. Using the symmetry 9 — > n — 9 we then arrive 
to an approximation to the ratio 5 n (x; M) = ^Mjg) given by 



^Erfc 2 (ixVIOf)4 A, 1 + t 2 V 1 " M - 1 "l 2sm »l 

At the next step we assume that the ^—integral is dominated by the finite values 
< 9 < 7r/2. This allows to replace the lower limit of integrations by and also to 
expand he ~ 1 + ln |2sin#|/lnM for InM — > oo. Changing then the integration variable 
r = (1 + w/lnM) 1//2 and keeping only the leading order terms the u— integral is easily 
calculated and the emerging integral takes the form (4 sin 2 9) x ^ d9. The latter 
is convergent as long as x < a/2 where it can be reduced to the Euler beta-function, 
see ( !T4|) . As a result, we arrive to the following expression: 

r n — r 2 /2) 

J im Snte M ) = T^Ti — hr\ ~ < x < ^ ( 29 ) 

Af^OO 1 Z [1 — X Z /4:) 

which is fully equivalent to the variance result ( 122|) we have anticipated on the basis of 
the thermodynamic formalism. For x > a/2 the 9— integral diverges at the lower limit 
9 — > rendering our large-M procedure invalid. This case will be separately treated in 
the end of the section. 

Now we are in a position to check numerically the range of applicability of the 
approximations derived. First we attempt to compare the results of exact numerical 
evaluation of the discrete sum in ( [271) to the integral ( 1281) . Actually, the direct evaluation 
of the sum in Mathematica is affordable up to M = 50000. To go up to the higher 
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Figure 5. Numerical simulations (red bars) compared to the theoretical prediction 
(127)1 (black line) for M = 64 (top) and M = 4096 (bottom). 1// signal is generated 
via the Fast Fourier Transform (FFT) method. 



values of M we use the identity 

M-x „ x f-i ,./, :i _ j 



/i -i rn„+i 

f{r)dr=Y,mn) / f(r)dr 
ir- i n = l 



(30) 



where R{n) = n for n < M/6 and R(n) = n — M/3 + | for n > M/6. Since for large 
InM the difference \h n — h n+ i\ is very small, the integral in the above expression can 
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be very accurately approximated by the trapezoidal rule 

2 / f(r) air « (h n+l - h n ) (f(h n+l ) - f(h n )) . 
This trick allows us to evaluate the sum in (ETJ) numerically up to M ~ 10 
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Figure 6. Comparison of the convergence to the asymptotic values. The partition 
function variance (blue circles , S z (q = 0.5, M)) reaches the asymptotic value 
S z (q = 0.5, oo) = 0.180341 for M ~ 2 16 . The relative variance S n (x = 1,M) for the 
counting function is presented as evaluated by three different methods: Filled Squares 
stand for the result of precise integration of the discrete sum ([27]) up to M = 50000. 
Empty Squares stand for a numerical integration involving the identity (|30[) and the 
trapezoidal rule. Diamonds describe the continuum approximation Eq. (|28p . The 
asymptotic value predicted by ([29]) for M — > oo is again 6 n (x — l;oo) = 0.180341 is 
nowhere close to the finite- M results and the agreement even seems to worsen with 
growing M. See the next figure. 



The results are presented figure 5 and figure 6 for x — 1 . In figure 5 we compare the 
fast convergence of S z (q = 0.5, M) against the very slow convergence of 5 n (x = 1; M) 
: for M ~ 2 15 5 z (q = 0.5, M) is already very close to the asymptotic value whereas 
8 n (x = 1 , M) shows a curious non-monotonicity and even for M ~ 2 30 we are still very 
far from the asymptotic value predicted by (129|) which is 5^\x = 1) = 0.180341. We 
observe that for M larger than 2 17 the continuum approximation of Eq.( |28l) matches 
perfectly the numerical integration involving the identity fl30|) and the trapezoidal rule. 
The latter matches perfectly the exact discrete sum (1271) even for moderate M. As we 
now are confident in the accuracy of the continuum approximation formula ( 1281) given 
by a double integral with InM entering as a simple parameter we can use it to check 
what values of InM actually ensure the validity of the asymptotic (1291) . The results are 
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presented in figure 6 and show that one needs astronomically big values of M even to 
achieve a rather modest agreement with the asymptotic value. 
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Figure 7. S n (x = 1,M) given by the continuum approximation formula (f28l) . 
The non-monotonic region seen in the previous figure is reflected in a tiny dip 
confined to a region close to the origin. The asymptotic value predicted by ([29| 
is S n (x = l;oo) = 0.180341. 

These facts explains our failure to confirm the infinite-M asymptotic by direct 
simulation of the variance of the counting function predicted for the periodic 1/f noise 
in the framework of the thermodynamic formalism. In any realistic signal analysis of 
such variance one must therefore rely upon the exact formula ( J27l) or its analogues 
instead of the asymptotic value ( 1291) . The conclusion should be of significant practical 
importance, in particular in view of the growing interest in numerical investigations of 
statistical properties of high values of the modulus of the Riemann zeta-function and 
of the characteristic polynomials of large random matrices. 

Having verified by the independent approach the asymptotic of the second moment 
of the counting function for x < \/2 we can now apply a similar methods beyond that 
range. The formal divergence of the right-hand side in ( 129]) for x — > \[2 simply means 
that the second moment of the counting function is not proportional to the squared 
typical scale Af e , but is parametrically larger and is actually determined by the detail 
of the probability density ( !l~9l) close to its upper cut-off n c which is not yet known to 
us. An accurate analysis of the second moment in the range a/2 < x < 2\/2 performed 
in the Appendix B shows that the leading asymptotic behaviour of S n (x; M) is given 
by 



5 n {x-M)^x\l ^^ M r M^- l )\ v / 2<x<2v / 2 (31) 
[2 - X/V2J 
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As is easy to see 5 n (x; M) \^Nm(%)j in the above domain, and is therefore the 
dominant term in the relative variance ( 125]) . The above result can be used to estimate 
the upper cutoff parameter n c . Indeed, as 5 n (x; M) ^> 1 the formula ()25l) in the range 
< x < 2 can be rewritten n 2 « V (\ — 5 n (x; M), with 5 n (x; M) given by ( 13~T1) . 

2-4 

On the other hand, in that range of x the probability density ffl9|) implies n 2 ~ n c 1 . 
Comparing the two results and using f[3"Tj) yields the scaling of the upper cutoff: 

n c ~M^^, v / 2<a;<2 (32) 



5. The position of threshold of extreme values in generic 
disorder-generated multifractal patterns 

Results obtained so far in the paper suggest a natural question about statistics of 
high values and positions of extremes of more general power-law correlated multifractal 
random field with a general non-parabolic singularity spectrum. Most obvious examples 
include the variety of the Anderson transitions[14j, but in fact many more random 
critical systems should be of that sort|10j. A straightforwards calculation outlined 
in [22] shows that behind each pattern of such type lurks a certain logarithmically 
correlated field, though in general of a non-Gaussian nature. Below we sketch that 
simple argument for the sake of completeness. Consider a d— dimensional sample of 
linear size L, and assume following [10] that the multifractal patterns of intensities p(r) 
is self-similar 

P q (r 1 )p s (r 2 ) oc L~ yq - S \r 1 - r 2 \~ Zq '% q,s>0, a < |r x - r 2 | < L (33) 
and spatially homogeneous 

p<j( ri ) = _L J where L = [ p q (r) dr oc L~ Tq (34) 

L J\r\<L 

The consistency of the above two conditions for |ri — r 2 | ~ a and |ri — r 2 | ~ L implies: 

Uq,s — d + T q+S , Zg tS = d + T q + T s — T q+S (35) 

If we now introduce the field V(r) = lnp(r) — lnp(r) and combine the identity 
■^p s \ s =o = lnp with the fact that r = — d we straightforwardly arrive at the relation 

V(r 1 )V(r 2 ) = -g 2 In ^ - r ^ , g 2 = -^-T q+S \ s=q=0 (36) 

Thus we conclude that the logarithm of any multifractal intensity is a log-correlated 
random field. The above argument does not say anything about higher cumulants 
of the field V(r), but it is easily checked that had those fields be always Gaussian 
the resulting singularity spectrum f(a) obtained from r q via the Legendre transform 
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would be invariably parabolic. Therefore, any non-parabolicity of the singularity 
spectra necessarily implies non-Gaussian nature of the underlying log-correlated fields. 
Nevertheless, combining our previous insights with results for disorder-generated 
multifractal patterns in [H] suggests the way in which our results on Gaussian 1/f 
noise can be generalized to statistics of high values and positions of extremes of more 
general non-Gaussian logarithmically correlated random processes and fields. 

As was already mentioned in the introduction, in the case of the Anderson 
transition the probability density of the inverse participation ratios I q depends only 
on the scaling ratio z = I q /I q \ with with I q standing for the typical value. Moreover, 
that ratio is expected to be power-law distributed: V g (z) ~ z~ x ~^ q [HIE]. We may try 
to combine that fact with the theory developed in the present paper to conjecture the 
typical position of the extreme values (maxima or minima) in a pattern of multifractal 
probability weights pi ~ M~ ai i = 1, . . . M ~ L d . A brief account of such a procedure 
is as follows. Suppose the mean participation ratios are given by I q = B(q) M~ Tq , 
with a coefficient B(q) of order of unity, and concentrate on those q for which typical 
and annealed exponents coincide. From it we recover in the usual way the singularity 
spectrum f(a) by the corresponding Legendre transform: f(a) > for a G [a_, a+] and 
further assume a_ > to avoid complications related to the so-called "multifractality 
freezing" [9j [131 12H] which would require a special care. Define a(q) as the solution of 
equation q = f'(a) and denote z q = J °° V q (z) z dz. Further given any function <f) q of the 
variable q define a "Lagrange conjugate" function 0*(ai), by the relation 0* («(<?)) = 4> q . 
Then we suggest that the density of exponent — Sf=i ^ (iff/ — a ) should be 

given asymptotically by the following "improved" multifractal ansatz: 



where n = n*(a) is a random coefficient of the order of unity distributed for a 
given a according to the probability density V*(n) defined in terms of V q {z) via 
the rule V*^(n) = V g (n). Indeed, substituting such an ansatz into the definition 
= Jo°° M~ qa pM {oi) da and performing the integral by the Laplace method for 



is distributed according to the probability density V q (z). This is precisely what is 
required, provided we identify I q ^ = I q /z^. Now we can substitute the same ansatz to 
the definition of the counting function N < (a) = J_° ra /)(a) da choosing a to the left of 
the maximum of f(a). This gives asymptotically N < (a) = n*(a)N t (a) where the scale 
Nt(a) is now given by 




(37) 




where the random variable z 



n*{a{q)) 



N t (a) 



z*(a)f'(a)\ 27rlnM 



B,(a) /|/"(«)| 



a_ < a < «o 



(38) 
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and defines the typical value of the counting function. A few typical "extreme" values 
among p^s in the multifractal pattern will be of the order of p m = M~ am , where a m is 
determined from the condition that the counting function becomes of the order of unity: 
Nt{oi m ) ~ 1. Clearly, at the leading order a m = a_ given by the left root of f(a) = 0, 
and the goal is to extract the subleading term. For doing this a crucial observation 
taken from [14J is that for q — >■ q c = /'(a_) the tail exponent u q of the IPR probability 
density V g (z) ~ z ~ 1 ~ UJq should tend to u qc = 1. As the derivative ou' q=gc is generically 
neither zero nor infinity the average ~Zq~ diverges as Yq~ ~ (q c — q)' 1 . In turn, as z*(a) 
is the Lagrange conjugate of 1q~ the divergence of the latter implies similar behaviour 
z*(ce) ~ (a — a_) _1 in the vicinity of a_. At the same time generically f'(a),f"(a) 
neither vanish nor diverge at a = a_ , and we do not see any reasons to expect that -B* (a) 
vanishes or diverges at this point either. Approximating f(a m ) ~ /'(a_)(a m — a J) we 
arrive at the following equation for the extreme value threshold a m : 

f'(a-)(a m - a-)hxM- -lnlnM + In (a m - a J) =0 (39) 
solving which for In M ^> 1 to the first non-trivial order gives 

a m ^a^ + U ^-\npj^ ^ a „lnM+^— ^lnlnM (40) 

2 j («_) hiM 2 J (a-) 

which is our main prediction for the typical position of the threshold of extreme values 
in disorder- induced multifractals. In particular, the value of the absolute maximum 
Pmax will be such that y = \np max — In pm is a random variable of the order of unity. 

6. Discussion and Conclusion 

In conclusion, we have studied, both analytically and numerically, the strongly 
fluctuating multifractal pattern in structure of those values of the ideal 1/ / nose which 
exceed a given high level comparable with the absolute maximum of the signal. The 
exploitation of the thermodynamic formalism allowed us to translate the distribution 
of the partition function found in the previous studies [301 ED] to a similar distribution 
for the counting functions. The latter is then characterized by a power-law forward tail 
giving rise to a parametric difference between the mean and the typical value of the 
counting function when the level position approaches threshold x m of extreme values. 
This difference, which can be traced back to the logarithmic correlations inherent in 
the 1/f noise explains the universal coefficient of the subleading term in the position 
of the threshold x m . 

We have also performed a direct numerical simulations of the 1/f signal and 
calculated numerically the lowest two moments of the partition function. This served 
to demonstrate that for the samples of M ~ 10 6 points the numerics follows M = oo 
results rather faithfully. Performing the same check for the counting function moments 
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however showed that truly asymptotic results can be never achieved even with very 
moderate accuracy due to prohibitively slow convergence. Instead, even for M as big 
as M ~ 10 9 one has still use a more elaborate finite-M formulas which we derived 
in the present paper for that goal. This lesson may prove important in view of the 
growing interest in numerical simulations of related systems arising in the framework 
of the Random Matrix Theory and the Riemann zeta function along the critical line 



Finally, by comparing the results obtained for 1/f noises in our paper with those 
known to hold for multifractal patterns of wave-function intensity at the points of 
Anderson transitions [9], Ej we propose a quite general formula for the position of 
extreme values in generic disorder-generated multifractal patterns with non-parabolic 
singularity spectra. We hope that such prediction can be checked against the accurate 
numerical data in random multifractals of various origin, and will generate further 
interest in statistics of high and extreme values in such multifractal patterns. We leave 
this issue as well as a related, but much more difficult question about actual statistics 
of the counting function in the region of extreme states for future investigations. 
For 1/f noise the latter should involve understanding of how the so-called "freezing 
phenomena" known to have profound influence on the behaviour of the partition 
function Z q with \q\ > 1 are reflected in the thermodynamic formalism correspondence 
between Z q and the counting function. Note that the freezing mechanism suggested in 
[291 EH1 EE] predicts for \q\ > 1 the tail behaviour for the distribution of Z = Z q to 
be V\ ql>1 {Z) ~ Z-( 1+ i)lnZ. It is based on considering properties of the generating 
function g q (y) = exp{—e~ qy Z q /Z e (q)} which in the limit M — > oo is conjectured to 
stay q— independent (i.e. "frozen") to the value g q =i(y) for all q > q c — 1. The factor 
In Z in "P|g|>i(Z) plays a prominent role and is believed to be a universal feature within 
the class of Gaussian logarithmically correlated fields. Note that such factor will be 
precisely absent for i.i.d. case. To that end it is curious to observe that in the context of 
the Anderson Localization [HI |9] a certain use of the thermodynamic formalism for the 
IPR's combined with a clever heuristic power counting lead to the distribution for I = I q 

(i ■ q c \ 
...... j ...... , ,, ,, , , , i J. It is most natural to suspect 

that the logarithmic factor In I may be present in the V q>qc (I) as well, and the accuracy 

of the procedure used in [HI [9] was simply not enough to account for it. Closely related 



questions are whether the generating function g q {y) = exp|— e qy I q /I q j will be 
actually q— independent for q > q c = f'{aJ) for general non-parabolic multifractals and 
whether the probability density of the logarithm of the (appropriately shifted) absolute 
maximum y = In pm P ' — lnp max , with prn given by (HOl) . will show a characteristic 
non-Gumbel tail |2/|e~' y ', y — > — oo [291 EB EI] as our extended analogy would suggest. 
All these intriguing issues certainly deserve further investigation, both numerically and 
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Appendix A: Variance of the counting function for Gaussian sequences. 



Suppose we have a sequence of M correlated Gaussian-distributed variables 
Vi,...,Vm characterized by the common variances V- 2 = Co and covariances VjVj = 
c iji i 7^ 3- Define 

Am (a) = J2 - a ) d ( V i " a ) ~ Wi " «) ' Wj " «)] (-1) 
Our goal is to show the validity of the following expression 

A M(a) = iexp(-f)f f'j^e-^', hv-M? (.2) 
For proving it we need the following 

Proposition. Suppose Vi,V 2 are two Gaussian-distributed variables characterized by 
the common variances V\ 2 = V± = Co and the covariance ViV 2 = c. Then for any two 
functions fi(V) and /2(V) with finite means fi, 2 (v) holds the identity 

d 



dc 



' fi(Vi) f 2 (V 2 j\ = f[(V l )f^(V 2 ) (.3) 



where /' = jj^. To verify the proposition we introduce the vector v = ^ j^ 1 
C = I C ° C ] and write the joint probability density of V\ and V 2 as 



denote 



c c 



p( v ) = e vT ° lv = r r c -u T c4>+i(v^ 1+ v 2( i> 2 ) #i#2 ^ 

2W detC J-™J~™ 2n 



from which it is immediately clear that 
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This implies: 

§- c [fi(V)f 2 (V)] = j°° j°° hiV,) h {V 2 )^ v V{w)dV l dV 2 = WoW) (.6) 

where the last equality follows after integration by parts. 

Now applying the proposition to the particular case /i(V) = f 2 (V) = 9(V — a) 
gives in view of jtw9(V — a) = 5{V — a) the relation 



d_ 

dc 



1 



9(V 1 -a)9(V 2 -a) = V{V l = a,V 2 = a) = — -e"^ (.7) 

2ix^4 - c 2 



and since 9(V% — a) = 9(V 2 — a) is independent of c the integration shows that 



D 12 = 9{V X - a) 9(V 2 -a)- 0(V X - a) ■ 9(V 2 - a) = - / —?== e ~ 

n Jo a/ c\ - u 2 



C An, „2 



2 r du 



Finally, introducing the variable r = y we convert the above expression to the 
form 

4 / a 2 \ f 1 dr _i£ r 2 /c - c , . 



which after applying to each pair of V*, Vj in (fTTJ) immediately implies ( T2"j) . 
Appendix B: Large-M asymptotic of the 5^ (x) for < x < 2. 



Our goal is to extract the large-M asymptotic behaviour of the integral featuring 
in OIL that is: 



For x > a/2 we anticipate that the main contribution for M — > oo comes from 9 — > 
and rescale the integration variable as 9 = nM~ u , u e (0, 1). Such a rescaling allows us 
to write h(9) — >■ = « / so that we have 



JO J hi', 



d r -In Mr 2 - 2 
1 + r 2 



e ~ In M a; 2 T 2 /4 / 2^ 

h(u) 



= TfM" 1 ^L_ e -lnM,^V4 _1 f 1 duh > {u){l + u) e -lnM(« + 4^) 

io i + r 2 y 

For InM ^> 1 the first integral is dominated by the lower limit and yields the 
leading order contribution JmL>i( x ) ~ TohVTnM- Second integral is dominated by 
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the vicinity of the minimum of the function C(u) = u + achieved at «* = -7= — 1. 

For a/2 < x < 2a/2 we have w* G (0, 1) so applying the Laplace method. Using 
= — 1 + \/2x — x 2 /A and -^C(u = u*) = 2\/2/x we find the leading order 
contribution given by: 



r 3 - -' 

2 v / 21nM(2 v / 2 



- J tt^ .1 S . M 1 -^ (.3) 



Finally, using asymptotic formula Erfc f|\/hi Mj ~ x 2 ^. g=ig we arrive at the expression 
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